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Abstract 



The nuclear binding energies for 28 nuclei including several isotopic chains with masses 
ranging from ^ = 64 to A = 226 were evaluated using the Skyrme effective nucleon-nucleon 
I interaction and the Extended Thomas-Fermi approximation. The neutron and proton 

■ density distributions are assumed in the form of Fermi functions the parameters of which 
. are determined so as to minimize the total binding energy of any given nucleus. The 

present study is restricted to quadrupole shapes, but the neutron and proton density 

■ distributions are free to have different deformations. A simple expression for the variation 
. of the nuclear energy with the neutron-proton deformation difference is derived. 

^ . PACS number:21.60.Jz,21.10.Dr,21.10.-k,21.10.Pc 

> 

^ : 1 Introduction 



The analysis of the experimental data on electron and a-particle scattering, pionic atoms, 
and annihilation of antiprotons show that, in most nuclei, neutrons and protons have slightly 
different r.m.s. radii The main reasons for this difference are the Coulomb repulsion between 
protons and unequal numbers of neutrons and protons. 

Fully microscopic approaches of the Hartree-Fock (HF) type using the effective nucleon- 
nucleon interactions of the Skyrme or Gogny type f^, ^, as well as the relativistic mean field 
theory reproduce in a rather satisfactory way the experimental proton and neutron r.m.s. 
radii and their isotopic shifts @, |^, |^. 

In deformed nuclei, neutron and proton density distributions are expected to have not only 
different radii, but also exhibit different shapes, i.e. different quadrupole and higher multipole 
deformations. An analysis of theoretical densities based on a surface multipole moment expan- 
sion shows that significant differences between neutron and proton deformations often occur 
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both in the ground state as well as along the whole path to fission p. The Hartree-Fock- 
Bogolubov (HFB) calculation performed with the Gogny effective force in Ref. for ^^^Th, 
238-u g^^^ 240p^ have shown that the multipole deformations of the proton and neutron 
density distributions of fissioning nuclei are far from being equal. The relative difference be- 
tween them exceeds often 10 % and undergoes large variations which means that the thickness 
of the neutron skin does not remain constant as the fissioning nucleus elongates. The effect on 
the nuclear binding energy of these deformation differences was found to reach approximately 
1.5 MeV, with fluctuations of the order of 1 MeV 0. These variations are clearly not negligible 
compared to typical flssion barrier heights. 

Until now a large majority of macroscopic-microscopic type calculations of potential energy 
surfaces of fissioning nuclei assume equal deformations of proton and neutron distributions for 
both densities and single-particle potentials (see e.g. p], ^ [l^l)- In view of the results obtained 
in Ref. , one can expect that such calculations predict fission barriers that are systematically 
about 1 MeV too high. In a sens things are even worth in the case of macroscopic-microscopic 
calculations as compared to the selfconsistent Hartree-Fock approach where deformation en- 
ergy surfaces are usually generated through calculations with a constraint on some multipole 
moments of the mass distribution, which leaves protons and neutrons free to deform differently 
within this constraint, even though it of course turns out that these deformations are quite close 
to one another. One must note, however, that only the fluctuating part of the correction for 
the proton and neutron deformation difference can actually influence barrier heights calculated 
with the macroscopic-microscopic approach, since the average value of this difference can be 
taken into account in the fltting procedure of the parameters of the macroscopic (e.g. liquid 
drop) model. Neglecting the neutron-proton deformation differences in both the macroscopic 
and microscopic parts of the potential energy could then lead to an overestimation of sponta- 
neous flssion lifetimes by a few orders of magnitude in heavy and superheavy nuclei and could 
also have a non negligible effect on the binding energy difference between neighboring nuclei 
what will be reflected in predicted lifetimes for the a or electron-capture decays. 

The aim of the present investigation is to develop a new term in the liquid drop type mass 
formula which will approximate the average variation of the binding energy when protons and 
neutrons deform in a different way. It seems to us that the most suitable approach to this 
goal is the Extended Thomas-Fermi (ETF) approximation |TT] in connection with the Skyrme 
energy functional . After a short outline of the model we present results of our calculations 
for a large sample of nuclei and derive a simple approximate expression which nicely reproduces 
the effect of the proton and neutron deformation difference in a deformed macroscopic model. 



2 Description of the model 

2.1 Skyrme density functional and the Thomas— Fermi approxima- 
tion 

In order to determine nuclear binding energies one often uses effective nucleon-nucleon inter- 
actions of the Skyrme type which describes quite accurately nuclear ground-state properties as 



well as the low-lying collective excitations |]T2[. For such an interaction the total energy density 
£{r) is an algebraic function of the neutron and proton densities pn and Pp, of the kinetic energy 
densities and Tp and the spin-orbit densities J„ and Jp 
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+ + t2)pr+ i(t2 - h){Pn Tn + PpTp) + ^(^2 " 3ti)p (1) 

+ ^(3tl + t2)(Pn VVn + Pp^'^Pp) + ^hPPnPp 



J- Vp+Jn- Vpn + Jp ■ Vpp] + ^(^1 - t2){Jp + C^'r 



where non-indexed quantities such as p stand for the sum of the corresponding quantities for 
neutrons and protons, as e.g. 

P = Pn + Pp (2) 

and where to, ^i, ^2, ^3, xq and Wq are force parameters which can be adjusted so as to reproduce 
nuclear ground-state properties (binding energies, radii, nuclear spectra, etc.). 

In a first approximation the nuclear energy (without Coulomb term) could be calculated by 
using the Skyrme functional (0) and the Thomas-Fermi (TF) approximation where 

r, = l{Mf''pT , q = {n^p] . (3) 

There is no contribution at this (TF) order of the semiclassical expansion to the spin-orbit 
density Jq as the spin has no classical analogon. These terms are therefore omitted in the 
following calculations. 



2.2 The model with sharp-surface distribution 

If only ellipsoidal deformations are considered, one can describe the surfaces of nuclei by a very 
simple parameterization, using a single deformation parameter a. The lengths of the axis of 
the axially symmetric ellipsoides are then given by 



a = b = i?Oq exp 
c = i?o, exp [crj 



1^ 

2^^ 



{n,p} (4) 



where and Rq*^ are the radii of the neutron and proton distributions respectively. Obviously 
a > correspond to prolate and cr < to oblate deformations. Introducing the isospin 
parameter 

the following dependence on A and / was obtained in the RMF calculation of Ref. 

i?o, = 1.237 (^1 - 0.157/ - ^) A'/' fm , (6) 

i?o„ = 1-176 fl + 0.250 /+ ^4^) A^/=^fm. (7) 



A 

For ellipsoidal deformations it can be shown that the nuclear surface is given by 

Rg{e) = . (8) 

/a2 cos2(6') + c2 sin2(^) 
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It is obvious that in the case of different deformations and sharp distributions of protons 
and neutrons there will appear regions in which there exist only protons or neutrons and other 
regions where both types of particles coexist. In order to evaluate the nuclear energy in such 
a case, the Skyrme functional depending on the nucleon densities, can be separated in three 
terms depending only on the proton density, only on the neutron density and an interaction 
term depending on the densities of both type of particles : 



1 3 

-{tl + ?,t2)pq Tq + ^ (^2 " tl)pq Pg 



(9) 



where q — {n,p}, 



Vpn = to{l + Y)ppPn + ^{tl+t2){ppTn + pnrp) 

- ^(3^1 - t2)(PpVVn + Pn^^Pp) + \h{plpp + pIPu) 



Finally, 8 can be written as 



(10) 



The total nuclear energy Enuc{N, Z) is then given as the volume integral oi S{r). 

Enuc{N,Z)= jE{f)dh . (12) 

In the case of a sharp surface distribution the changes of the nuclear part of the energy as a 
function of deformation can be just expressed as the product of the volume Vtpn where protons 
and neutrons coexist and the term of the energy density connected with the interaction of these 
two types of particles : 



where 



pn 



^-^nuc{^(^^ ^pn^pn j 



1 - 



exp(3A(j) — 1 

:g^)^exp(-Aa)- 
1 — exp(— SAcr) 



3/2 



- 1 



3/2 



(13) 



(14) 



with Act = (T„ — o". 

One should notice that for this type of nucleon densities, the terms of S{r) with Ap vanish, 
so that the nuclear energy doesn't change during the deformation process as long as the proton 
and neutron distribution don't cross each other, since the volumes between the sharp surfaces 
is kept constant. After the surfaces intersect the energy rapidly decrease because there the 
region in which protons and neutrons coexist Qpn changes its volume as a function of the de- 
formation. This fact leads to the changes of the term Vpn- while terms Vp and Vn are constant. 
A typical behavior of the binding energy with growing proton-neutron deformation difference 
is presented in Fig. 1. The model with the uniform density distribution is, of course, too rough 
to make a realistic estimates but nevertheless it should give some idea about the main effect. 
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2.3 The model with the Fermi-type distribution. 



More accurate estimates of the change of the binding energy when the protons and neutrons 
deform differently can be made with diffuse density profiles. In the following we have chosen 
for our analysis the density distributions in the form of Fermi functions 

where po, are the saturation density parameters obtained from the normalization conditions : 

Pp(r)rfV = Z , J pn{r)dh = N (16) 

and where Rq{0) is given by Eq. (^. The surface width parameters a^, are 9 dependent and 
equal to 



ag{e) = ^sm\e) + cos^e) + sm\e) cos2(e)[exp(3(7g) + exp(-3(Tg)] ao, (17) 

a definition which guarantees that the surface diffuseness is constant along the direction per- 
pendicular to the surfaces of the ellipsoids. The parameters {i?o„, ao„i -Rop, OOp} characterizing 
the spherical density distribution of neutrons and protons are obtained by a minimization pro- 
cedure of the energy functional (0) within the second order ETF approximation [|ll| and using 



the Skyrme SIII force . 



3 Results of the fits 



The calculation were performed for 28 even-even nuclei involving several isotopic chains from 
Ni to Th. We suppose that the deformation of these nuclei can be characterized by a single 
global deformation parameter a plus a quadrupole type deformation parameter A/3 = (3n — $p 
measuring the difference in the proton versus neutron deformation. The global deformation of 
the nucleon distribution can be e.g. defined through the parameter a introduced by Myers and 

g 

n 2 



Swiatecki 



a 



27r / de 







R{e) - n 



00 



00 



where R{0) is the half-density radius of the nucleon distribution of the deformed nucleus and 
i?oo the radius of the corresponding spherical distribution. For an ellipsoidal deformation as 
characterized by Eq. @) a can be easily expressed in terms of the parameter a and one obtains 
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a 



(19) 



In order to be able to express the total energy as a function of the difference A/5 between the 
proton versus neutron deformation we need to find a way to express this difference in a way 
which is independent of the specific choice of the deformation parameters like e.g. in Eq. @). 
We have chosen to define parameters I3q as : 



V20 

Voo 



(20) 
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where the Q^^q are the monopole and Q2Q the quadrupole moments of the proton and neutron 
distribution with : 

^?S = J P.if^r'Pxd'r , (21) 

and where 

Voo ~ Voo "T Voo 5 

Qoo = Qoo* ■ ^ ; (22) 

Voo ~ Voo ■ ^ • 

are the so-called average monopole moments for protons and neutrons respectively. 

We have found that contrary to the rapid change (Fig. 1) of the binding energy in case 
of the sharp-edge density distribution the nuclear part of energy varies almost parabolically 
with the proton and neutron deformation difference A(3 when the smooth density profile (p!5|) 
is assumed for protons and neutrons. 

The binding energy can then be parametrized in the following form 

B{A, I, a, Pn, h) = BavM, /, a) + a(A, /, a) ■ (^„ - 4)^ , (23) 

where Bam is the part of the energy that is generated for equal deformations for protons and 
neutrons and which can be obtained e.g. form the standard liquid drop or droplet model. To 
determine the stiffness parameter a of the parabolic form in A/5 we proceed in the following 
way : 

For various nuclei with different values of the mass parameter A and isospin parameter / and 
for a 7^ one finds that this parameter a{A, /, a) can be expressed in the form 

a(A,J,a) = ao(A^) [1 -C2«(l + C3/)] , (24) 

where the deformation independent coefficient ao(A, /) can be parametrized as 

ao{AJ) = cr[A^'^Y . (25) 

The values found by our fitting procedure for the parameters ci, C2, C3 and n are the following 

ci = 1.21 MeV , C2 = 3.0 , C3 = 3.0 , n = 4 . 

The stiffness parameter a{A, /, a) evaluated for all considered isotopic chains is compared on 
the l.h.s. of Fig. 2 at a = deformation with the parameter ao(v4, /), Eq. ([25|) . The a 
deformation dependence of the parameter a{A, J, a) is illustrated on the r.h.s. of Fig. 2. We 
have also compared the above estimate of 0(^4, /, a) with the results obtained using two other 
Skyrme interactions. We have found that the estimates obtained with the Skyrme SII and SVII 
forces are identical to those of presented above for SIII within a 10%. 

The quality of the fit of the binding energy as function of the deformation difference is 
demonstrated on Fig. 3 for three nuclei (^^Zr, ^^^Nd and ^°®Pb) and for three different values 
of the global deformation parameter a and compared with the corresponding ETF results. 

In order to test the predictive power of equation (^3]) we have performed an additional 
calculation for three dysprosium isotopes (neutron deficient, /3-stable and neutron rich) and 
compare in Fig. 4 their ETF energies (crosses) with the approximation given by ( ^3|) (solid 
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line) and this for three different deformation. As can be seen in Fig. 4 the agreement of both 
results is rather satisfactory in all cases, except for the very deformed neutron deficient ^^''Dy 
isotope. 

When a 7^ and {/3n — (3p) = 0, the energy is reproduced by the standard liquid drop model 
or any other macroscopic model. 



4 Discussion and conclusions 

We have found a very simple procedure to describe the influence on the total nuclear energy 
arising from the fact that neutron and protons might deform differently. The additional term 
in the nuclear mass formula which we obtain is proportional to the square of the deformation 
difference of of both distributions and increases as 

This suggest that it might be useful to generalize the currently used macroscopic-microscopic 
approaches in order to allow protons and neutrons to have different multipole deformations. 
One could start from a Strutinsky type prescription : 

EstUiP^^^}, {Pt^}) = E^UiP?}, + }) + SEl^UiPP}) , (26) 

where both the macroscopic and and microscopic components depend on proton and neutron 
deformations. Then, performing a minimization with respect to e.g., the neutron multipole 
deformations while keeping proton deformations as independent variables (or the other way 
around), should give a more realistic description of potential energy surfaces. Details of such a 
method as, for instance, the way to generalize macroscopic models in order to include different 
deformations for protons and neutrons, are left for a future study. 

It was already shown in Ref. ||^ that taking into account the proton versus neutron defor- 
mation difference could shift the position of the ground state and saddle point in the potential 
energy surface by as much as 0.5 MeV. Such an effect could improve the predictions of nuclear 
masses as given by macroscopic-microscopic type of models and change the predicted barrier 
heights. In consequence the effect discussed in the present paper could change the theoretical 



predictions of a-decay and spontaneous-fission life times H, 10]. 
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Figure captions 



1. Typical behavior of the binding energy with growing neutron-proton deformation differ- 
ence for the nucleus ^osp]^ g^j^^ ^-j^g sharp proton and neutron density distributions. 

2. The mass number (l.h.s. figure) and the deformation (r.h.s. figure) dependence (crosses) 
of the stiffness parameter 0(^4, /, a) (Eq. ^Bp) and its approximation by the formula (^) 
(solid lines). 

3. The change of the binding energies due to different neutron-proton deformations as ob- 
tained in the ETF approach is compared with the approximate expression of Eq. (^). 

4. The approximation (p3|) of the binding energy (solid line) made for three Dy isotopes and 
at three deformation points is compared with the ETF results (crosses) as function of the 
neutron-proton deformation difference. 
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